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1.  Introduction 

For  certain  models  of  viscoelastic  fluids  with  fading  memory,  classical  steady  channel 
flow  does  not  exist  beyond  a  maximal  wall  shear  stress.  This  occurs  when  the  shear  stress 
for  steady  flow  decreases  with  strain.  For  generalized  Newtonian  models  of  viscoelastic  flow, 
such  a  decrease  implies  that  the  flow  is  unstable.  Because  of  this  example,  a  constitutive 
relation  that  exhibits  a  maximal  wall  shear  stress  is  regarded  as  defective. 

-We'  report  on  work  showing  that,  contrary  to  this  intuition,  such  models  correctly 
describe  the  experimentally  observed  "spurt”" phenomenon:  exceeding  this  critical  stress 
results  in  a  large  increase  in  volumetric  flow  rate.  The  transition  to  spurt  flow  is  analogous 
to  a  dynamically  generated  phase  transition. 

To  analyze  this  phenomenon,  we  derive  a  system  of  conservation  laws  that  govern  the 
flow;  these  equations  take  the  form  of  gas  dynamics  with  relaxation  terms.  We  solve  the 
Riemann  problem  for  this  non-strictly-hvperbolic  system,  and  incorporate  this  solution 
into  the  random  choice  method.  Numerical  simulation  of  channel  flow,  with  the  maximal 
wall  shear  stress  exceeded,  shows  that  a  discontinuity  forms  at  the  wall,  allowing  the  fluid 
to  “slip”;  no  steady  state  exists.  However,  when  a  small  Newtonian  viscosity  is  included 
in  the  model,  a  “slip  layer”  forms  and  the  flow  approaches  a  discontinuous  steady  state. 
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2.  Viscoelastic  Fluids  with  Fading  Memory 

The  motion  of  a  fluid  is  defined  by  the  spatial  position  x(X,t )  at  time  t  of  each  fluid 
particle  X.  The  velocity  of  a  particle  is  then  v  =  dx/dt.  The  deformation  of  the  fluid 
is  characterized  by  the  deformation  gradient  F  =  dx/dX.  Alternatively,  the  deformation 
undergone  by  the  fluid  from  time  s  to  time  t  may  be  specified  using  the  relative  deformation 
gradient  F3(t)  =  F(t)F(s)~ 1  =  dx(t)/dx(s).  We  will  work  in  the  Eulerian  picture,  so  all 
fluid  variables  are  regarded  as  functions  of  the  spatial  position  x  at  time  t. 

When  the  flow  is  incompressible  and  isothermal,  the  equations  governing  the  evolution 
of  F  and  v  are  the  continuity  equation, 

F  =  (Vv)F ,  (2.1) 

and  the  conservation  of  momentum, 

p{T  =  V  •  a  .  (2.2) 

Here  p  is  the  mass  density,  a  is  the  Cauchy  stress,  and  a  dot  denotes  the  convective  time 
derivative;  also,  the  components  of  V£T  are  (Vv)*j  =  dv'/dx* .  We  will  write  a  =  —pi  +  E, 
where  p  is  an  isotropic  pressure  that  maintains  incompressibility;  E  is  the  extra  stress. 

The  response  characteristics  of  the  fluid  are  embodied  in  the  constitutive  relation 
for  the  stress.  For  viscoelastic  fluids  with  fading  memory,  these  relations  specify  the 
extra  stress  as  a  functional  of  the  deformation  history  of  the  fluid.  Many  sophisticated 
constitutive  models  have  been  devised;  see  Ref.  [1]  for  a  survey. 

In  the  present  work  we  focus  on  the  Johnson-Segalman  model  [5]  as  a  prototype  for 
general  constitutive  models.  This  model  accounts  for  non-affine  deformation  of  Gaussian 
networks  by  introducing  a  parameter  a,  —  1  <  a  <  1,  that  enters  a  nonlinear  generalization 
of  the  classical  Maxwell  model: 

E  =  f  m{t  —  s)  {Fs(t)Fs(t)T  —  f}  dfi  ,  (2.3) 
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where  Fs  satisfies 

P,  =  -RlF,  ,  F,(s)  =  I.  (2.4) 

Here  Ra  =  uj  -  ad,  with  d  =  4  [V£T+  (Vu)r]  and  u  =  §  [Vt?  -  (Ve)T]  being  the  defor- 
mation  rate  and  vorticity  tensors.  When  a  =  l,  Fs  =  Fs  and  the  Johnson- Segalman  model 
reduces  to  the  Maxwell  model. 

The  function  m  is  a  memory  kernel,  defined  on  [0,oo),  that  we  assume  to  be  smooth, 
nonnegative,  and  integrable.  A  typical  example  is 

m(s)  =  n\e~Xa  ,  (2.5) 

where  /i  is  a  shear  modulus  and  A  is  a  relaxation  rate.  With  this  choice  of  memory  kernel, 
the  Johnson- Segalman  model  may  be  formulated  as  the  differential  equation 

E-f-E/2a  +  i?rS  =  2a/id-AS  (2.6) 

for  the  extra  stress. 

Constitutive  relations  such  as  Eq.  2.6  exhibit  a  mixture  of  elastic  and  viscous  behavior. 
This  may  be  seen  heuristically  as  follows.  In  the  short  relaxation- time  limit,  A  — >  0, 
Eq.  2.6  shows  that  an  objective  time  derivative  of  £  is  proportional  to  the  deformation 
rate:  £  ~  2 a/id.  This  is  characteristic  of  elastic  behavior,  and  leads  to  the  interpretation 
of  aft  as  a  shear  modulus.  By  contrast,  when  A,/i  — ►  oo  with  /i/A  fixed,  £  ~  2(afj./X)d ; 
thus,  the  model  displays  viscous  behavior  with  a/.i/ A  being  the  Newtonian  shear  viscosity 
coefficient. 

As  shown  in  Ref.  [5],  the  constitutive  relation  for  the  stress  derives  from  the  energy 

£  =  ^-f  m{t  —  s)Tr  [F3(t)F3(t)T  —  I)  ds  .  (2.7) 

2 aP  y-oo 

The  energy  satisfies  the  conservation  principle 

p(±v-v  +  £)’  =  V-(v-<r)-pr,  (2.S) 
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where 


r  = 


_2op  J  m'(*  _*)Tr  {F,(t)F,(t)T  -  1}  ds 


(2.9) 


represents  the  rate  of  dissipation  of  energy. 


3.  Shear  Flow 

Interesting  properties  of  viscoelastic  constitutive  relations  are  exhibited  in  simple  pla¬ 
nar  shear  flow.  With  the  flow  aligned  along  the  y- axis,  the  deformation  gradient  F  takes 
the  form 

F=(!?)’  (31) 
where  the  strain  e  depends  only  on  the  spatial  coordinate  x  and  the  time  t.  In  particular, 
the  flow  is  incompressible;  we  assume  p  to  be  constant  throughout  the  fluid.  Let  v  denote 
the  y-component  of  v  and  let  a  denote  the  shear  stress  T,zy  =  T,yz .  The  variables  v,  cr, 
Err,  and  T,yy  also  depend  only  on  x  and  t,  while  the  pressure  p  takes  the  form  p(x,  y,  t )  = 
Po(x,t)  —  f{t)y.  The  quantity  /  corresponds  to  a  pressure  gradient  driving  the  flow. 

In  this  notation,  the  continuity  equation,  Eq.  2.1,  reduces  to 


-  vr  =  0  , 


(3.2) 


and  the  y-component  of  conservation  of  momentum,  Eq.  2.2,  becomes 


pvt  -  crx  =  f 


%3) 


Furthermore,  the  x-component  of  the  conservation  of  momentum  equation  determines  p0 
through  the  equation  px  =  (EIX)X. 

The  stress  components  are  expressed  in  terms  of  the  history  of  the  strain  through  the 
constitutive  relation,  Eq.  2.3.  As  is  easily  determined, 


a  =  — 


m(t  —  s)sin(a[e(f)  —  e(s)])ds  , 


(3.4) 
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where  a  =  y/l  -  a2.  Similarly,  Eri  =  -£(1  -  a)A/*i  and  Eyy  =  j(l  +  a)Aflt  where 
Ar\  =  Eyy  —  EZI  is  the  principal  normal  stress  difference: 

2a  1 1 

-  s)  {1  -  cos(a[e(f)  -  e(s)])}  ds  .  (3.5) 

®  •/  —  oo 

Also,  the  energy  is  given  by  p£  =  i-A/j.  Notice  that  M\  >  0. 

In  the  special  case  when  m  is  given  by  Eq.  2.5,  Eqs.  3.2-5  imply  the  differential 
equations 


pvt  -  ax  =  f  , 

(3.6) 

crt  —  [a/j.  —  o2p£']  vz  =  —Xa  , 

(3.7) 

p£t  —  crvz  =  —A p£  . 

(3.8) 

Several  mathematical  results  are  known  for  this  quasilinear  system; 

see  Ref.  [10]  for  dis- 

cussion  and  references. 

•  The  system  is  strictly  hyperbolic  if  and  only  if 

p£  <  > 

a.1 

(3.9) 

and  the  change  of  type  has  been  associated  to  certain  instabilities. 

•  For  smooth,  hyperbolic  initial  data  with  sufficiently  small  variation,  a  classical  solution 
exists  for  all  time. 

•  Nevertheless,  smooth  initial  data  exist  for  which  the  solution  develops  singularities  within 
a  finite  time.  As  the  critical  time  is  approached,  the  spatial  derivatives  of  v  and  a  become 
unbounded. 

Thus,  the  fading  memory  acts  as  a  weak  dissipative  mechanism:  the  source  terms  in  the 
equations  serve  to  counteract  the  formation  of  singularities  from  sufficiently  smooth  data. 
When  singularities  do  form,  however,  Eqs.  3.6-8  are  no  longer  valid  because  the  products 
of  distributions  £vz  and  avz  are  ill-defined.  Instead,  the  physical  conservation  equations, 
interpreted  in  the  sense  of  distributions,  govern  the  flow. 
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4.  Equivalent  Conservation  Form 

Eqs.  3.2-3  and  3.4  express  physical  conservation  principles  that  remain  valid  past 
“blow-up.”  They  are  not  standard  conservation  laws,  however,  because  the  stress  is  given 
by  a  convolution  integral.  To  solve  these  equations,  we  first  derive  an  equivalent  system 
of  equations  that  is  more  amenable  to  analysis. 

This  system  is  obtained  by  introducing  z  >  0  and  8  through  the  equations 

2  cos(a$)  =  m  *  cos(ofe)  ,  (4-1) 

z  sin(a$)  =  m  *  sin(ae)  ,  (4-2) 


where  we  use  the  convolution  notation 


*  <t> 0)  =  J 


m(t  —  s)<f>(e(s))ds  . 


Notice  that  z  and  8  are  continuous  in  t  for  almost  every  x,  assuming  that  e  is  of  bounded 
variation.  Therefore,  in  calculating  the  distributional  time  derivatives  of  2  and  8,  we  may 
legitimately  multiply  by  functions  of  z  and  8.  Simple  algebra  shows  that 

zt  —  m( 0)  cos  (a[e  —  #])  +  cos(a8)  m  *  cos(ae)  +  sin(a#)  m!  *  sin(ae)  ,  (4.4) 

a  2  8t  =  m(0)  sin  (a[e  —  #])  —  sin(a$)  m!  *  cos(ae)  +  cos(a<?)  m!  *  sin(ae)  .  (4.5) 

Furthermore,  the  stress  components  are  given  as  functions  of  e,  2,  and  8: 

a  =  aa~1  z  sin  (a[e  —  8))  ,  (4-6) 

A'!  =  2aa-2  {/itot  -  2cos(a[e  -  0])}  ,  (4.7) 


with  fitot  =  f0°°  m($)ds. 

In  this  way,  the  constitutive  relation  3.4  is  transformed  into  Eq.  4.6  together  with  the 
conservation  laws  4.4-5.  The  convolution  integral  that  appeared  in  the  definition  of  a  has 
been  removed;  in  its  place  is  a  pair  of  evolution  equations  for  the  “internal  variables”  2 
and  8.  Because  the  stress  depends  on  e  and  8  only  through  the  combination  r  =  e  —  8,  6 
may  be  regarded  as  the  unrecoverable,  or  “viscoplastic,”  part  of  the  strain,  while  r  is  the 
“elastic”  strain. 
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5.  Nonlinear  Waves 


Let  us  now  restrict  to  the  case  of  a  single  relaxation  mode:  m(s)  =  pX  exp(-Xs).  If 
we  eliminate  6  in  favor  of  r  =  e  -  9  by  combining  Eq.  4.5  with  Eq.  3.2,  then  Eq.  3.2  for  e 
decouples  from  the  system.  Therefore,  we  obtain  the  conservation  laws 

rt  —  vx  =  -pXa~l  2 —  1  sin(or)  ,  (5.1) 

pvt  -  ct(t,z)x  =  /  ,  (5.2) 

zt  =  —X  {z  —  ficos(ar)}  (5.3) 

together  with  the  constitutive  relation 

a (r,z)  =  aa-1  2rsin(ar)  .  (5-4) 

In  the  present  isothermal  approximation,  energy  is  not  conserved  in  discontinuous 
flows.  Rather,  the  energy  equation,  Eq.  2.8,  serves  as  an  “entropy”  condition  for  weak 
solutions: 

P  [y^2  +  £(ri2)]*  -  <vf-  V£(r,z)  (5.5) 

where 

p£(r,  z)  =  aa~2  {p.  —  z  cos(ar)}  .  (5.6) 

Except  for  the  presence  of  source  terms,  Eqs.  5.1-3  are  identical  to  those  for  La- 
grangian.  homentropic  gas  dynamics,  with  r/p  identified  as  specific  volume,  —a  as  pres¬ 
sure,  and  z  as  a  function  of  the  entropy.  Corresponding  to  the  fundamental  thermodynamic 
identity  is 

pdS  =  a  dr  -f  pT dz  ,  (5-7) 

where  T  =  — p~laa~ 2  cos(or)  plays  the  role  of  “temperature.” 

However,  the  equation  of  state,  Eq.  5.6  or  Eq.  5.4,  would  be  rather  peculiar  for  a  gas; 
cf.  Ref.  [8].  The  physical  state  of  the  fluid  is  unchanged  if  r  is  replaced  by  r  +  2ir/a,  so 
state  space  is  cylindrical;  for  convenience,  we  require  that  — 7r/2  <  ar  <  37r/2,  along  with 
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2  >  0.  The  wave  speeds  are  ±c  and  0,  where  pc 2  =  da  I  dr  =  az  cos(ar)  defines  the  sound 
speed  c.  Therefore  the  system  is  strictly-hyperbolic  if  — tt/2  <  or  <  tt/2.  whereas  there  is 
a  complex-conjugate  pair  of  purely  imaginary  wave  speeds  when  tt/2  <  ar  <  3~/2.  The 
wave  mode  with  zero  speed  is  linearly  degenerate,  and  the  sound  modes  with  speed  ±c  are 
genuinely  nonlinear  in  the  hyperbolic  region  except  at  r  =  0. 

Despite  the  peculiarities  of  the  system,  we  have  determined  completely  the  structure 
of  scale-invariant  nonlinear  waves.  Such  a  wave  consists  of  a  sequence  of  elementary  scale- 
invariant  waves,  either  centered  discontinuities  or  rarefaction  waves,  connecting  constant 
states  on  the  left  and  right.  Discontinuities  are  required  to  satisfy  Liu’s  generalization  of 
Oleinik’s  entropy  condition,  which  guarantees  that  the  energy  dissipation  inequality  5.5 
holds.  This  admissibility  condition  is  equivalent  to  requiring  shock  waves  to  have  viscous 
profiles:  with  Eq.  5.2  augmented  by  the  Newtonian  viscosity  term  ( r]vz)z ,  admissible  shock 
waves  arise  as  limits  as  77  — >  0  of  travelling  wave  solutions. 

The  wave  structure  is  conveniently  depicted  with  a  wave  curve,  the  locus  of  states 
U  =  (r,  v,  z)  on  the  right  for  a  fixed  state  Ui  =  (17,,  vi,zl)  on  the  left.  Only  2  may  change 
across  waves  with  zero  speed,  so  these  wave  curves  are  trivial.  On  the  other  hand,  2  remains 
constant  across  waves  corresponding  to  the  characteristic  families  with  speeds  ±c,  so  we 
suppress  the  2  component  of  U  when  drawing  wave  curves.  Fig.  1  shows  the  wave  curve  of 
the  — c  family  for  a  representative  initial  state  U l\  the  wave  curve  of  the  opposite  family 
is  obtained  by  reflection  through  the  vertical  line  v  —  vi.  The  figure  was  produced  using 
a  computer  program  that  constructs  the  wave  curves  for  general  systems  of  conservation 
laws  [4].  Of  particular  interest  is  the  manner  in  which  the  wave  curve  corresponding  to  the 
—  c  characteristic  family  continues  past  the  boundary  of  the  elliptic  region  with  waves  that 
contain  rarefaction  and  shock  waves  of  the  +c  family;  this  is  possible  because  the  wave 
speeds  coincide  at  the  elliptic  region. 

With  the  wave  curves  known,  Riemann  initial-value  problems  may  be  solved.  We  have 
written  a  computer  program  that  solves  Riemann  problems,  and  have  incorporated  it  into 


the  Glimm-Chorin  random  choice  method.  We  refer  to  Ref.  [9]  for  a  detailed  discussion  of 
this  method.  The  Riemann  problem  solver  can  also  be  used  in  enhanced  computational 
methods  for  two-dimensional  flow  problems,  such  as  the  front  tracking  method  [2]. 
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6.  Steady  Channel  Flow 

As  an  application  of  the  results  just  described,  we  have  investigated  a  flow  of  a 
Johnson-Segalman  viscoelastic  fluid  in  a  channel  [7].  In  particular,  we  have  studied  the 
approach  to  steady  flow,  since  this  corresponds  to  rheological  experiments. 

Let  the  vertical  walls  of  the  channel  be  located  at  x  =  ±X2W'  The  flow  is  symmetric 
about  the  centerline  (x  =  0),  so  we  focus  on  the  left  half  of  the  channel,  —  <  x  <  0. 

Symmetry  requires  that  vx  =  0  at  the  centerline,  and  the  no-slip  boundary  condition 
imposes  that  v  =  0  at  the  wall. 

Steady  shear  flow  occurs  when  the  velocity  v  is  independent  of  time.  In  this  case,  the 
strain  rate  k  —  et  =  vx  is  constant  in  time,  so  that  the  stress  a  coincides  with 

CL  00 

^steady  («)  =  —  /  m(s)  sin(oAcs)  ds  .  (6.1) 

Q  Jo 

The  steady  solution  may  be  determined  by  integrating,  in  turn,  the  ordinary  differential 
equations  —  <7steady(*0i  =  /  and  vx  =  «,  with  the  boundary  conditions  that  k  =  0  at  x  =  0 
and  v  —  0  at  x  =  ~\w. 

Steady  solutions  may  fail  to  exist  if  crgteadyC^)  does  not  vary  monotonically  with  k. 
For  example,  when  m(s)  =  ^Aexp(  — As), 

a/z  olk/ A 

^steady(K)-  q  X  +  (uK/A)2  (6’2) 

attains  a  maximum  crmax  =  ^a^/a  when  qk/ A  =  1.  Classical  steady  flow  in  a  channel 
of  width  w  cannot  exist  if  the  pressure  gradient  /  exceeds  a  critical  value  /cr ;t  given  by 
ffmax  =  jU’/crit-  Notice  that  this  phenomenon  is  not  associated  with  loss  of  hyperbolicity, 
since  pS  —  |a/z/a2  at  the  stress  maximum,  well  within  the  bound  3.9. 

Traditionally,  the  non-existence  of  classical  steady  solutions  caused  by  a  shear  stress 
maximum  is  regarded  as  a  defect  of  the  constitutive  law.  This  conclusion  is  based  on 
intuition  appropriate  for  generalized  Newtonian  models  of  viscoelastic  fluids.  Shear  flow 
for  such  a  fluid  is  governed  by  the  single  equation 

pvt  -  [r](vz)vz]z  =  f  ,  (6.3) 
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corresponding  to  having  a  viscosity  coefficient  rj  that  depends  on  strain-rate.  In  a  flow 
regime  where  tj(k)k  decreases  with  strain-rate  k,  however,  Eq.  6.3  has  the  character  of 
a  backward  heat  equation,  which  suffers  from  the  Hadamard  instability.  Therefore,  for 
generalized  Newtonian  fluids,  t](k)k  must  increase  with  k  in  a  physically  stable  steady 
solution. 

The  Johnson- Segalman  model  has  the  same  steady  solutions  as  does  a  generalized 
Newtonian  fluid  with  t7(k)k  =  <7steady (*),  so  one  might  think  that  it  exhibits  the  same 
instability  in  regions  where  crstea(jy  decreases.  This  conclusion  is  not  warranted,  however, 
because  the  system  of  equations  maintains  its  evolutionary  character  even  near  the  shear 
stress  maximum. 

A  question  remains:  what  is  the  solution  if  the  critical  pressure  gradient  fcrit  is 
exceeded?  To  answer  this  we  must  solve  the  time-dependent  equations,  allowing  for  the 
possible  formation  of  singularities. 

7.  Slip  and  Spurt  Phenomena 

We  have  simulated  numerically  a  simple  experiment  for  a  fluid  governed  by  the  sys¬ 
tem  5.1-4.  Without  loss  of  generality,  units  are  chosen  so  that  a/z  =  1,  A  =  1,  and  p  —  1, 
and  the  fluid  variables  are  scaled  so  that  a  disappears  from  the  equations.  Given  this 
choice,  the  channel  thickness  w  and  the  pressure  gradient  /  (which  depends  on  time)  axe 
the  essential  parameters  defining  the  flow.  In  the  numerical  experiment,  the  channel  width 
is  w  =  1,  so  that  the  wall  of  the  channel  is  located  at  x  =  —  |  while  the  center  is  at  x  —  0. 
The  flow  is  initially  in  the  classical  steady  state  corresponding  to  the  critical  pressure 
gradient  /cr ;t;  then  the  pressure  gradient  is  increased  abruptly  to  the  super-critical  value 

1.2/cnf 

The  result  of  this  calculation  using  the  random  choice  method  is  shown  in  Fig.  2.  The 
fluid  velocity  v  is  plotted  vs.  position  x  at  successive  time  intervals;  generally  the  velocity 
increases  with  time.  During  the  early  stages  of  the  experiment,  the  flow  settles  into  a 
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quasi-steady  state.  This  effect  is  especially  evident  in  a  plot  of  the  centerline  velocity  as 
a  function  of  time,  and  it  is  more  pronounced  when  the  channel  width  w  is  smaller  [7]. 
Eventually,  however,  a  thin  layer  develops  at  the  wall  in  which  the  velocity  rises  to  a  value 
that  is  nearly  constant  across  the  channel.  For  practical  purposes,  the  fluid  has  broken 
free  from  the  wall  and  is  accelerating  uniformly  under  the  applied  pressure  gradient;  thus 
the  fluid  “slips.” 


X 

Fig.  2:  Onset  of  slip  for  a  fluid  without  Newtonian  viscosity. 

To  understand  this  result  better,  the  same  experiment  was  performed  for  a  fluid  that 
has  a  small,  but  nonzero,  Newtonian  viscosity  (e.g.,  because  of  the  presence  of  a  solvent). 
For  such  a  fluid,  the  term  (r; vx)z  is  added  to  the  right-hand-side  of  Eq.  5.2.  So  long  as 
0  <  t]  <  1/9  (in  dimensionless  units),  the  steady  shear  stress,  plotted  vs.  strain,  still 


exhibits  a  maximum. 


Fig.  3  shows  the  results  for  tj  =  0.01,  as  calculated  using  the  Lax-Wendroff  method 
with  Tyler  artificial  viscosity.  Now  the  slip  layer  is  much  thicker,  with  its  outer  boundary 
corresponding  to  a  discontinuity  in  the  strain  rate  vx.  The  solution  approaches  a  steady 
state  in  which  vx  is  discontinuous  but  the  total  stress  Steady  (ur)  4-  t]vx  is  continuous. 


-0.5  -0.4  -0.3  -0.2  -0.1  0 

X 

Fig.  3:  Onset  of  spurt  for  a  fluid  with  Newtonian  viscosity. 


Further  study  of  channel  flow  for  a  Johnson-Segalman  fluid  yields  several  observa¬ 
tions  [7]. 

•  The  nominal  shear  rate  (based  on  volumetric  flow  rate),  plotted  against  the  nominal 
shear  stress  (based  on  the  pressure  gradient),  shows  an  abrupt  increase  at  a  critical  point, 
i.e.,  “spurt.” 

•  The  thickness  of  the  slip  layer  (or,  equivalently,  the  stress  at  which  the  strain-rate 
discontinuity  occurs)  exhibits  hysteresis  as  a  function  of  the  loading. 
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•  There  are  three  essential  time  scales  to  the  flow  problem,  corresponding  to  processes  gov¬ 
erned  by  Newtonian  viscosity,  by  wave  propagation,  and  by  relaxation.  In  some  situations 
this  leads  to  a  quasi-steady  supercritical  flow  that  persists  for  a  time  orders  of  magnitude 
longer  than  the  wave  propagation  scale. 

There  are  two  essential  features  responsible  for  spurt  and  hysteresis  phenomena  in 
viscoelastic  models:  (1)  a  non-monotone  stress-strain  relation  in  steady  flow;  and  (2)  a 
dynamical  mechanism  for  selecting  steady  shock  waves  that  are  physical.  For  example,  the 
spurt  and  hysteresis  phenomena  were  analyzed  by  Hunter  and  Slemrod  [3]  for  a  related,  but 
simpler,  model  of  a  viscoelastic  fluid.  Their  analysis  shows  how  the  first  two  observations 
may  be  understood,  and  we  expect  that  it  may  be  applied  to  the  present  model.  We 
emphasize,  however,  that  the  Johnson-Segalman  model  differs  in  an  essential  way  from  the 
model  considered  by  Hunter  and  Slemrod:  it  remains  hyperbolic  even  in  the  range  where 
the  stress  decreases  with  strain,  as  is  desirable  for  equations  that  describe  evolution  of  a 
physical  system. 

Spurt  and  hysteresis  offers  an  explanation  of  certain  observations  in  rheological  exper¬ 
iment  [11].  Quantitative  agreement  between  experiment  [11]  and  Johnson-Segalman  model 
has  been  verified  in  Ref.  [6].  Therefore,  we  have  confidence  that  the  interesting  mathemat¬ 
ical  structure  of  models  of  viscoelastic  flow  reflects  observable  physical  phenomena. 
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